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Abstract 



In order to analyse the lattice dependence of ferromagnetism in the 
two-dimensional Hubbard model we investigate the instability of the 
fully polarized ferromagnetic ground state (Nagaoka state) on the tri- 
angular, honeycomb and kagome lattices. We mainly focus on the local 
instability, applying single spin flip variational wave functions which 
include majority spin correlation effects. The question of global insta- 
bility and phase separation is addressed in the framework of Hartree- 
Fock theory. We find a strong tendency towards Nagaoka ferromag- 
netism on the non-bipartite lattices (triangular, kagome) for more than 
half filling. For the triangular lattice we find the Nagaoka state to be 
unstable above a critical density of ra = 1.887 at U = oo, thereby sig- 
nificantly improving former variational results. For the kagome lattice 
the region where ferromagnetism prevails in the phase diagram widely 
exceeds the flat band regime. Our results even allow the stability of 
the Nagaoka state in a small region below half fllling. In the case of 
the bipartite honeycomb lattice several disconnected regions are left 
for a possible Nagaoka ground state. 

Keywords: Hubbard model; Ferromagnetism; Triangular lattice; 
Kagome lattice. 



1 Introduction 



The Hubbard model is well established today as the simplest model to de- 
scribe strongly correlated electrons on a lattice. It contains hopping processes 
for each spin direction (one orbital per lattice site) and an on-site repulsion 
as a simple approximation to the Coulomb interaction between the electrons. 
One of the most important motivations for introducing the Hubbard model 
thirty years ago [1, 2, 3] was to obtain a better understanding of the appear- 
ance of ferromagnetism in transition metals. However, the investigation of 
the existence of ferromagnetic phases in the Hubbard model, in particular 
the proof of the stabihty of the fully polarized ferromagnetic state, the so- 
called Nagaoka state, turned out to be a very hard task. Up to now only very 
few exact results are known about the existence of a Nagaoka ground state 
in the Hubbard model. Most of them are extensions of the famous Nagaoka 
theorem [4] and deal with the very special case of a half-filled lattice plus 
one single electron at infinite on-site repulsion U (for example [5, 6]). In this 
situation the Nagaoka state is found to be stable on most of the common 
lattices like square lattice, sc-, fee-, &cc -lattice or triangular lattice. For bi- 
partite lattices this stability holds also for the particle-hole symmetric case 
of one hole in an otherwise half-filled band. 

The first thermodynamically relevant proof of the stability of the Nagaoka 
state was presented by Mielke and Tasaki [7, 8, 9] for a special class of 
two-dimensional lattices with fiat energy bands, for example the kagome 
lattice. For all other lattices progress was only achieved in restricting the 
region where the Nagaoka state might be a ground state. Starting with the 
work of Shastry, Krishnamurthy and Anderson (SKA) [10] the instability of 
the Nagaoka state was investigated using single spin fiip variational wave 
functions [11, 12, 13]. This work focussed mainly on the square lattice. Two 
of the present authors achieved an improvement of former results by including 
majority spin correlation effects into the variational ansatz [14], again for the 
square lattice. 

As demonstrated by Mattis [15] two criteria are essential for the Nagaoka 
mechanism to produce ferromagnetism: there should be loops in the under- 
lying lattice (which is not fulfilled for the one-dimensional chain e.g.) and 
there should be constructive interference of different paths of particles or 
holes around loops. Like other non-bipartite lattices the triangular lattice 
shows the shortest possible loop length three and is therefore a good candi- 
date for Nagaoka ferromagnetism in the case of more than half filling. For 
less than half filling, however, the Nagaoka state is expected to be suppressed 
because the second criterion is violated. 
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In the first part of the present paper we use various methods to investigate 
the instabihty of the Nagaoka state on the triangular lattice. In order to 
obtain a general view of the phase diagram we start with the Hartree-Fock 
approximation (section 2.1). In this context we investigate ferromagnetic 
states as well as other spin structures like the 120° configuration with three 
sublattices. Allowing the possibility of phase separation we achieve a large 
reduction of the Nagaoka stability region near half filling. In section 2.2 a 
systematic analysis of the stability of the Nagaoka state with respect to a 
single spin flip is performed by applying some of the extended variational 
wave functions of [14] to the triangular lattice. This requires an eflRcient 
method for calculating the elements of the one particle density matrix, which 
we describe in detail in the Appendix. We are able to show that the Nagaoka 
state is not stable for all densities above half filling at C/ = oo. For one 
hole less than half filling the Nagaoka theorem is known to be not valid for 
non-bipartite lattices. We study this case in more detail in section 2.3 and 
demonstrate that a global Nagaoka instability follows from rather general 
considerations which apply not only to the triangular lattice. 

In sections 3 and 4 we continue our study of Nagaoka ferromagnetism on 
two-dimensional lattices with the honeycomb and the kagome lattice. The 
densities of states of both lattices are closely related to the DOS of the tri- 
angular lattice. Again we use variational single spin flip wave functions to 
establish regions of guaranteed instabihty of the Nagaoka state in the phase 
diagram. For the bipartite honeycomb lattice two types of regions are left for 
a possible Nagaoka ground state in which different physical mechanisms ap- 
ply. The non-bipartite kagome lattice has been the subject of growing atten- 
tion in recent years. The interest arises primarily because a two-dimensional 
kagome like structure is found in real substances [16]. Furthermore, the ex- 
istence of a flat energy band in the nearest neighbour hopping model gives 
rise to interesting physics. For example, it implies that the Nagaoka state 
is stable in the flat band region since a ferromagnetic alignment of the spins 
requires no additional cost of kinetic energy and guarantees the absence of 
double occupancies. However, as Mielke [7] pointed out there is a degener- 
acy of the ground state, including states with all possible values of the total 
spin. Nevertheless, using statistical arguments he showed that for a certain 
density range states with extensive total spin dominate independently of the 
strength of the Coulomb interaction U. In our work we go beyond the region 
considered by Mielke deriving variational criteria for the instability of the 
Nagaoka state in the whole range < n < 2 and for all possible values of U. 
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2 Triangular lattice 



By diagonalizing the hopping part of the conventional single band Hubbard 
model [17] 

n = -t + n^]n,l (1) 

on the triangular lattice we obtain the band dispersion 

£A(k) = -t (2 cos{K) + 4 cos j cos ky^ j (2) 

where k belongs to the hexagon-shaped first Brillouin zone (BZ). For t > 
the lower band edge £h = —z^t = — 6t is found at k = 0, whereas the 
upper band edge with the energy = St is reached at all the corners of the 
hexagon. The density of states Pa(£) = (<5(£ — £a{^)))bz (expressed by a 
complete elliptic integral in the Appendix) is depicted in Fig. 1. It shows 
a logarithmic van Hove singularity at e = 2t where the Fermi surface forms 
a hexagon with a volume of exactly 3/4 of the whole BZ and contains the 
saddle points of the band (2) (e.g. Q = (0, 27r/A/3) ) as corners. 

In the phase diagrams displayed in this paper we will always represent 
the on-site repulsion U in terms of U^ed = U/{U + Ubr) where Ubr denotes 
the Brinkman-Rice critical coupling [18]. For the triangular lattice Ubr = 
15.81 \t\. 



2.1 Phase diagram in HF-approximation 

In this subsection we present a phase diagram of the triangular lattice in 
HF-approximation, thereby extending previous studies of [19, 20, 21]. In 
particular, we obtain a first instability region for the Nagaoka state which 
will be extended in the following subsections. 

Let us first briefly review HF-theory: In (unrestricted) HF-theory the two 
particle interaction is replaced by 

Un,^n,^ ^ ({Sn,)n, - 4(S,) S, - i( {6n,y - 4(8,^ - n^)\ , 

i i 

where n stands for the averaged particle density, Sui := Ui — n denotes the 
charge density fluctuations and corresponds to the spin operators at site i. 
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This ansatz simplifies if one restricts the theory to a uniform charge den- 
sity {{6ni) = 0) and to a hehcal spin density wave (S^) = ^ (cos(QRj)ea, + 
sin(QRi)e^ ) with m denoting the local magnetization. Under this restric- 
tion Fourier transform and diagonalization lead to the variational energy per 
site 

EsDw = J E(^-(k) n_(k) + E4k) n+(k)) + ^ + ^r^^ (3) 

k 

where we have used the abbreviations 

E±(k) = e (k) ± Y^A2 + Ae2(k) 

£(k) ] g(k + Q/2)±g(k-Q/2) 
Ae(k) J 2 

and where we introduced A for Um/2. 
Minimizing the energy with respect to A results in 

U 2L^ ^A2 + A£2(k) 

which is equivalent to the fulfilment of the self-consistency condition m = 
2|(S)|. In the limit A ^ the right hand side of this equation represents 
the paramagnetic susceptibility x(Q) and the critical U at which the para- 
magnetic state becomes unstable against a second order transition to a helical 
configuration with twist vector Q is given by 

1 1 v^ n(k-Q/2)-n(k + Q/2) 

Ucr £(k + Q/2)-£(k-Q/2)' 

We now turn to the special case of the triangular lattice for which we 
obtain the phase diagram presented in Fig. 2. The boundary between the 
paramagnetic phase and helical configurations is indicated by the lower full 
line. We find that the instability of the paramagnetic state occurs for those 
Q-vectors for which the two Fermi surfaces e(k ± Q/2) = ep are just tan- 
gent to each other (imperfect nesting). With increasing particle density, the 
favoured Q-vectors follow straight lines in the BZ along the para-helical phase 
boundary as depicted in Fig. 3: starting at the center Q = of the BZ for 
£f = —Qt (point A in Figs. 2 and 3) the Q vector reaches Q = (0, 2Tr/\/S) on 
the boundary of the BZ for £p = —2t (point B), moves along the boundary 
and reaches the corner Q = (27r/3, 27r/\/3) of the BZ for £p = —t (point 
C) and then turns back towards the center as ep = 2 t is approached. In 
the limit of ^i? ^ 2 t (which corresponds to n ^ 1.5) the Stoner criterion 
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tells us that the critical coupling for Q = approaches zero due to the 
logarithmic van Hove singularity: U^^^ = Pa{£f) oc — 3/(47r^t)ln |2 — Sp/t]. 
However, this ferromagnetic instability is not the leading one, because the 
Fermi surface is a perfectly nested hexagon. The corresponding nesting vec- 
tors Q = (0, 27r/\/3), (±7r, tt/'s/S) all correspond to alternating rows of up 
and down spins. This structure will be called the Neel structure in what 
follows. 

Due to this nesting the favoured Q vector right before £p = 2 t jumps 
from a value near zero to one of the nesting vectors for which Ucr then van- 
ishes as tU~^ oc (In |2 — i.e. it drops faster than corresponding to the 
Stoner criterion. For sp > 2 t the favoured Q first moves from one of the 
nesting vectors towards Q = 0, jumps onto the line A-C for even higher £p 
and finally returns towards Q = 0. 

For Coulomb interactions above the instability of the paramagnetic state 
we find areas in the {U^n) phase diagram where the optimum Q varies con- 
tinuously as well as other areas where the optimum Q has a fixed value such 
that particular spin patterns are stable in a whole area. The stable patterns 
are 

i) the ferromagnetic spin pattern Q = 0, 

ii) the antiferromagnetic three sublattice 120° structure Qaf = (47r/3,0), 

iii) the Neel structure of alternating rows of up and down spins QNeei = 
(0,27r/V^). 

Let us now briefly discuss the difl'erent stabihty regions starting with the 
antiferromagnetic 120° structure. At half filling this structure is favoured 
for all U exceeding 5.27 t (see [20, 21] for a more detailed discussion). In 
the limit of large U the stability of this structure follows easily from the 
asymptotic expression for the SDW energy 

+ £(Q)t Zt'^ + g(QAF)t _ Wt_ 

2C7 - W ~ ~ 2Z7' 

_ zt'^ + £(Q)t 

with W denoting the band width {W/^, = 9t). 

The 120° alignment remains favoured if the density is lowered moderately 
(n < 1). To see this we start at half filling and remove one single particle {U 
shall be large enough to avoid band overlap). The energy changes by 



EsDw{n = 1) = 
m = 
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dE = min(-£;_(k) n) = min(- £ (k) + Ja^ + A£^{k) ) 

k 2k' 2 

> -3t-|(l-m). (6) 

Equality only holds for the ferromagnetic structure and for the 120°-structure. 
As the latter structure has the lowest magnetization for given U it is not 
only favoured at n = 1 but also gains most energy if one particle is removed. 
Therefore the AF is stable in a finite range of densities n < 1 and the favoured 
twist vector doesn't change with n. Decreasing n further we find a first order 
transition to the Nagaoka state (for large U), whereas for lower U there is a 
first order transition to the Neel-structure. 

If, on the other hand, the density is increased above half filling the 120° 
structure immediately gets unstable and the favoured Q continuously turns 
towards Q = 0. These intermediate helical states are however not really 
stable for large U as the corresponding energy is not convex as a function of 
the density n. Therefore the system is unstable against a phase separation 
between the antiferromagnetic and the Nagaoka phase and one has to use the 
Maxwellian construction to restore convexity. The critical density Tier below 
which phase separation takes place in dependence of the Coulomb coupling 
follows from the equation 

EsDw{U, n = l) = Ej^{n„) - E'j^{n„){n„ - 1), (7) 

where Ej>^ denotes the energy of the Nagaoka phase. The asymptotic eval- 
uation for [/ ^ oo of (7) leads to tlU = Air {n^r - lf/{sVs). This is the 
strongest reduction of the stability range of the Nagaoka state in the limit 
U ^ oo, we were able to achieve and will not be exceeded in the following 
sections. Especially, it is clearly superior to the linear asymptotic behaviour 
t/U oc Tier — 1 which would be the result of analysing the instability of the 
Nagaoka state against a vanishing twist d'^Egow /dQ'^lQ^o = 0. This last 
instability gets relevant, however, if one lowers U as the corresponding curve 
crosses the phase separation curve discussed above. Below the crossing point 
phase separation takes place between the 120° structure and helical states 
Q = {qi,0) and the instability line of the Nagaoka state is determined by 
d^EsDw/dQ%^o = 0. 

Of course, the Nagaoka state might also be unstable against a partially po- 
larized ferromagnetic state. If depolarization is supposed to be a continuous 
process it can be studied in terms of a single spin flip problem. The obvious 
criterion is that the gain of band energy due to the flipped spin £t — £f should 
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equal the cost in Coulomb energy U{2 — n). The explicit calculation shows 
that this HF single spin flip instabihty is never the relevant criterion, i.e. 
we found no regions in the phase diagram where a less than fully polarized 
ferromagnetic state is stable. 

So much for the HF phase diagram of the Hubbard model on the tri- 
angular lattice as far as hehcal SDW states are considered. Summarizing 
our results concerning the Nagaoka state one has to realize that its stability 
region is clearly overestimated by the HF approximation. In particular, HF 
theory doesn't exclude stability for any filling factor in the limit U ^ oo. 

A well-known ansatz for improving HF theory uses a Gutzwiller projector 
controlling double occupancy. Unfortunately, except for high dimensions 
[24, 25] a rigorous evaluation of the energy of Gutzwiller wave functions is 
not possible. It is, though, possible for the very special but important case 
of a Nagaoka state with only one spin flipped. Gutzwiller single spin flip 
wave functions were used by SKA [10]. The results of this work will be the 
starting point of our next subsection. 

2.2 Variational single spin flip states 
2.2.1 Formulation 

The variational single spin flip states which will be discussed in the following 
were formulated in order to achieve an energy lower than the Nagaoka energy 
and thereby prove the instability of the Nagaoka state in a certain density 
range. The general method consists of minimizing the energy of the spin flip 
state \SF) with respect to the variational parameters which are contained in 
the trial wave function. From the condition 

(SF\n\SF) ^ 

with E// being the Nagaoka energy, one attains an area in the phase dia- 
gram where the Nagaoka state is deflnitely unstable. Whenever Ae > the 
Nagaoka state might be a ground state. 

According to the structure of the Hubbard Hamiltonian Ae contains three 
parts: the change in kinetic energy of the spin up particles due to the spin 
flip, the kinetic energy the down spin gains by delocalizing and the potential 
energy which results from double occupancies made possible by the spin flip. 
Under the restriction that the particles involved in the spin flip are given 
definite momenta the lowest possible energy is achieved by removing the 
spin up particle from the Fermi edge and creating the spin down particle 
with a momentum qj belonging to the lower band edge ej. 
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For our calculations we always assume n < 1. This is possible because the 
particle hole transformation maps the case n > 1 to less than half filling by 
changing the sign of t. In the phase diagrams for the non-bipartite lattices, 
however, we show the full density range < n < 2 for positive t. 

2.2.2 SKA Gutzwiller wave function 

The variational ansatz originally considered by Shastry, Krishnamurthy and 
Anderson [10] consists of a Gutzwiller projector applied to a single spin flip 
state constructed as explained above 

\SKA) = n(l-(l-5)«iT«u)<aCk,T lA/') (9) 
1 

= ;^E^^^''(i-(i-pKKck,TlA/') 

with the variational Gutzwiller parameter g controlling the probability of 
double occupancies. For g = 1 correlations are not taken into account and 
we recover the rather trivial single spin flip problem already mentioned in 
the HF part of our paper. The case U = oo, however, corresponds to ^ = 
when double occupancies are completely prohibited. 

The spin flip energy Aeoo for infinite U is calculated to be [10] 

Aeoo = e| + e| , e-^ = -{ej>//6) - ep , ei = £b 6 {1 - {ej>/ / 6 ztf) (10) 

with the lattice coordination number z, the Nagaoka energy per lattice site 
and the hole density 6 = 1 — n. There is no potential energy due to the 
absence of double occupancies. Evaluating Ae = for finite U by optimizing 
g for each given U and S leads to an analytic expression for the Nagaoka 
instability line Ucr{S) [14]: 

uUS) = - ''^'^'-l - '^^^^ . (11) 

(1 - djAeoo 

The lattice structure enters into (11) via the density of states which is needed 
to calculate 6 and e^/- as functions of the Fermi energy. 

In Fig. 4 we show the instability line obtained for the triangular lattice. 
As expected the Nagaoka state turns out to be completely unstable for less 
than half filling, whereas for more than half filling the spin flip energy remains 
always positive for large U . In order to demonstrate the importance of the 
Gutzwiller projector we also show again the HF result from Fig. 2 obtained 
by setting g = I which leads to a gross overestimate of the stability range of 
the Nagaoka state, especially for less than half fiUing. 
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2.2.3 Including majority spin hopping processes 

An extended single spin flip ansatz was introduced by Basile and Elser in 
1990 [11] to investigate the instability of the Nagaoka state at infinite U. 
In order to lower the kinetic energy loss of the spin up particles due to the 
existence of the flipped spin they extended the SKA Gutzwiller wave function 
including spin up hopping terms. These terms allow the spin up particles to 
move away from the down spin position 1 to other lattice sites r instead of 
simply projecting out double occupancies: 

|5^°') = ;^Ee'^''ciTE/r-lCrVt Ck.TlAT) . (12) 

The amplitudes of the hopping processes are controlled by the variational 
parameters fr i which are assumed to reflect the symmetry of the underlying 
lattice. 

Basile and Elser investigated this trial state numerically for a finite square 
lattice with 61 sites. Hanisch and MuUer-Hartmann [14] used a modification 
of this ansatz making it suitable for an evaluation in the thermodynamic 
limit and for all values of U. Starting from nearest neighbour hopping they 
successively added further spin up hopping processes. In order to consider 
finite U double occupancies were made possible by including the term g ■ 
with the Gutzwiller parameter already known from the SKA wave function 
(9). The Basile- Elser like ansatz for nearest neighbour hopping only is given 

by 



1 



(1 - (1 + /ici|5]c+^ 



cJck^tlA/') ■ (13) 



The amplitude of the hopping processes to the neighbouring sites n of the 
down spin is controlled by the variational parameter fn_i = /i, the varia- 
tional parameter /o has been set to 1. The lattice structure enters into the 
calculation of the spin fiip energy for \BEnn) again only via the density of 
states. As a consequence this ansatz allows the investigation of the instabil- 
ity of the Nagaoka state in the thermodynamic limit on an arbitrary lattice 
without great numerical effort. We will apply (13) also to the honeycomb 
and kagome lattices in sections 3 and 4, respectively. 

Going beyond this ansatz by including further hopping processes we ob- 
tain the spin fiip energy in general as a rational function in the elements 
of the single particle density matrix. In [14] Hanisch and Miiller-Hartmann 
evaluated these matrix elements for the square lattice. For the present cal- 
culations we have developed a method to calculate the elements of the single 
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particle density matrix for the triangular lattice efficiently such that only one 
numerical integration has to be performed for each matrix element. As de- 
scribed in the Appendix we used recursion formulae to express the integrand 
in terms of complete elliptic integrals of the first and second kind. The min- 
imization of the spin flip energy with respect to the variational parameters 
corresponds to the solution of a generalized eigenvalue problem. 

The long-dashed line in Fig. 4 shows the Nagaoka instability according to 
the Basile-Elser criterion for the triangular lattice including nearest neigh- 
bour hopping only. For more than half filling the SKA result is significantly 
improved, in particular in the region of large U and n. For the first time it 
becomes obvious that the Nagaoka state is not stable for all n > 1 at infinite 
U . The critical density above which the Nagaoka state cannot represent the 
ground state of the Hubbard model is calculated to be = 1.912. An even 
lower value of = 1.887 is achieved if further hopping processes are taken 
into account. This decrease of the critical density is achieved by including 
into the variational ansatz spin up hopping to the 6 first, 6 second and 12 
fourth neighbours. Including hopping to even more distant neighbours does 
not affect the result within the numerical accuracy as the importance of the 
processes diminishes with increasing distance from the fiipped spin. The 6 
third neighbours are an exception in this context since they do not infiuence 
the result for n^. 

As already mentioned the improvements achieved by the Basile-Elser 
ansatz in comparison to SKA are most pronounced near n = 2. In the t = —\ 
picture for less than half filling, which we used for our calculations, this 
corresponds to the range of large hole densities where the hopping processes 
are most eflRcient. For small hole densities most of the sites near the down 
spin position are occupied and therefore few hopping processes are allowed, 
especially for large U where this extension of the SKA ansatz therefore has 
only a small effect (see Fig. 4). In this region of the phase diagram, however, 
a distinct improvement of the SKA result is achieved by allowing the system 
to gain the antiferromagnetic exchange energy of the order jU between the 
down spin and the neighbouring up spins. This is done by replacing the 
Gutzwiller projector in (13) with another term allowing double occupancies, 
namely introducing a simple local correlation effect by binding a hole to the 
neighbours of the doubly occupied site: 

1 - (1 - g)n,^ (1 - ni^) + « E(l " "mt) • (14) 

m 

Here m stands for the nearest neighbour sites of the down spin site 1 and 
ct is a variational parameter. The full fine of Fig. 4 marks the Nagaoka 
instabihty line obtained with a variational single spin fiip state including the 
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Basile-Elser hopping processes mentioned before and the exchange term (14). 
As for the square lattice [14] this ansatz leads to a qualitatively better result 
than the Basile-Elser wave function with a simple Gutzwiller projector and 
correctly reproduces the divergence of the critical coupling Ucr in the limit 
n \ 1. The critical density n„ = 1.887 is not influenced by the inclusion of 
the exchange term while the minimal critical coupling is slightly improved to 
Ucr = 9.13 t which is the value at n = 1.586. 

We also considered another local correlation term introduced in [14] which 
makes the amplitude for the creation of a down spin at the lattice site 1 de- 
pend on the spin up occupancy of the nearest neighbours r of 1 in order to 
improve the spin down kinetic energy. This effect is included in the varia- 
tional ansatz by extending (14) 

(l-ni^) — (l-niT)(l+/5E(l-^rT)) (15) 

r 

with an additional variational parameter p. For the square lattice, this leads 
to a pronounced reduction of the critical hole density [14] compared with the 
SKA and BE results whereas for the triangular lattice it does not improve at 
all the value of Ucr obtained with the Basile-Elser hke ansatz. The physical 
reason for this is the fact that (again in the t = —1 picture for less than half 
filling) the hole densities being in consideration for the Nagaoka instability 
at infinite U for the triangular lattice are much larger than for the square 
lattice and consequently making an additional restriction by projecting out 
double occupancies at the neighbouring sites of the down spin position is 
inefficient since most of these sites are unoccupied anyway. For intermediate 
hole densities the term is also not relevant since the critical coupling U„ is 
so small that the system prefers the moderate potential energy caused by 
double occupancies rather than the loss of spin up kinetic energy due to the 
strict suppression of double occupancies near the flipped spin contained in 
(15). 

2.3 A static scattering problem 

In the limit of small hole densities the Nagaoka state was shown to be unstable 
for all flnite values of U on arbitrary lattices already in the early work of 
Richmond and Rickayzen [26]. In their variational single spin flip state the 
down spin is created at a fixed lattice site and this way represents a static 
scattering potential for the spin up particles. The kinetic energy of the 
fiipped spin is neglected which is a good approximation for large U and 
small 6 because under these conditions the spin down particle is unlikely to 
move. 



13 



The spin up particles no longer occupy Bloch states but scattering states 
and the calculation of the energy spectrum is reduced to a single particle 
problem. Richmond and Rickayzen used the Green's function method to 
derive an analytic formula for the energy difference between the Nagaoka 
state and their variational state: 

1 7 ( T^UpU) \ 

Ae = — de arctan — — \ — £p . 16 

7r4 \\-UA{£)j ^ ' 

Here A and vrp stand for the real and the imaginary part of the local prop- 
agator, respectively. A is calculated from the density of states p by using a 
Kramers-Kronig relation. 

Richmond and Rickayzen achieved their general results by evaluating (16) 
in the asymptotic limit oiU ^ oo. In our present work we determined the 
Nagaoka instability line for the triangular lattice according to the Richmond- 
Rickayzen criterion exactly by numerical evaluation of Ae = for all 1 < 
n < 2. We want to remind the reader that the limit n \ 1 is mapped to the 
Richmond- Rickayzen limit (5 \ by changing the sign of t. 

To evaluate (16) completely for all values of n a double numerical inte- 
gration is necessary. In purpose to diminish the numerical inaccuracies it is 
useful to treat singularities in an analytical way. As the density of states 
Pa(£) is known to diverge logarithmically at £ = 2t it is possible to find 
a regularized expression for pa- By inserting this regularized form of the 
integrand into the Kramers-Kronig relation one part of the integral can be 
expressed explicitly in terms of dilogarithm functions. This way we showed 
that the real part A of the local propagator has a discontinuity by at 
£ = 2t. The integrands of the remaining integrals which have to be solved 
numerically are all continuous functions. 

In Fig. 5 our result for the Richmond-Rickayzen state is compared with 
the results of other variational criteria which have turned out to be reasonable 
near half filling (i. e. n\ 1), namely the single spin fiip state (14) including 
exchange processes and the two curves obtained in HF theory by evaluating 
the instability of the Nagaoka state against an infinitesimal twist and against 
phase separation, respectively. The Richmond-Rickayzen criterion is the best 
local instability criterion in the limit 6^0 since on 6/\n{l/6). The 
global instability due to phase separation, however, leads to a substantially 
stronger criterion since it is the only one where the critical coupling diverges 
as U^^^ oc S'^. 

For larger hole densities the Richmond-Rickayzen ansatz is a poor vari- 
ational state due to the fact that there the kinetic energy of the spin down 
particle plays an important role in reducing the energy of the spin flip state. 
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Beyond 8{eF = 0) = 0.399 the Richmond-Rickayzen ansatz is not snfRcient 
for infering the instability of the Nagaoka state for any t/ > (see Fig. 5). 



2.4 Global instability of the Nagaoka state 

As we have already seen by evaluating the SKA ansatz the Nagaoka state is 
locally unstable for less than half filling on the triangular lattice. In partic- 
ular, for (1 — n) -C 1 and U = oo the SKA wave function leads to an energy 
gain of 

Ae = {l-n)[p-\e,)l2-zt{l-e'J{ztn 
= (l-n)t{7:V3-9)/2^-1.78t{l-n). 

In this section we will present a general consideration which shows that a 
single hole produces a global instability of the Nagaoka state. 

Let U = oo, t > and consider a lattice of L sites with periodic boundary 
conditions and N = L — 1 electrons (one hole). We are going to show that 
whenever there are at least two inequivalent points k in the BZ with e(k) = £t 
then there exist helical states with energy E < Ej^ = — e^. 

The finite energy gain dE = E — Ej^ (independent of L) makes sure that 
the instability also holds for (possibly very small but) finite hole densities 
in the thermodynamic limit. Therefore our proof is of thermodynamic rele- 
vance. As a by-product the explicit construction of such states provides in- 
formation about which magnetic structures are favoured by the hole hopping 
although we are not able to treat the paramagnet in an adequate manner. 

To prove the above statement we start from an arbitrary (HF) helical 
spin structure: 

IQ) = n (cos(^)c+ - sin(^)4) |0). 

To proceed it is useful to define holon and doublon operators in such a way 
that the above helical structure corresponds to a fermion vacuum: 



— )c+ - sm(— )c+, di = sin(^-)ci| + cos(^- 



The hopping part of the Hubbard Hamiltonian then reads (R^j := Rj — R^ 
Ht = -t E[cos(^^)(/i./i; + dfd,) + sm{^^){h,d, - dthi) + h.c.]. 
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Now one easily calculates the energy of the following superposition of one 
hole states 

|k>Q = h+\Q) = L-i/2^exp(ikR,)/i+|Q) 

i 

to be 

^^'^(k, Q) = -{e{k + Q/2) + e{k - Q/2))/2 = - ~e (k) 

which reproduces equation (6) in the hmit U ^ oo. Defining 
e^P(Q) = miiij^ '(k, Q) we obtain the energy s^i \o) = —£t for the Nagaoka 
state. If there are two band maxima at inequivalent points ki and k2 in 
the BZ, we obtain the same energy £^p{Q) = —£t for a helical state with 
Q = ki — k2 7^ 0. However, this state is not an eigenstate of the hole 
hopping Hamiltonian VTitV {V denotes the projector onto states without 
double occupancies). Applying the analytic Lanczos method to this state 
we can therefore construct helical one hole states with energies below —St- 
Performing the first Lanczos step explicitly we obtain the state 

I^)q = i-'/'Eexp(^kR.)/i+(a + E/5i-.<^^ )IQ) 

i i 

(flipped spin at site j neighboured to the hole at site i). Assuming the 
normalization laj^ + l/^i-^P = 1 we obtain the trial energy 

ef{k, Q) = S (k) - ^tsin(^) (exp(-M,,) a^l^ + c.c.) 

i 

QR 

j ^ p 

The sum over covers all common nearest neighbours of the neighboured 
pair of sites % and j. For the triangular lattice the number of such sites p for 
given % and j amounts to = 2. 

The evaluation of £^^'(k, Q) is especially transparent for the choice = 
/5exp(^(k — Q/2)Rij)/-y2, where one obtains 



e (k) 



A£(Q/2) \ 



„j Ag(Q/2) £(Q-k)+^i£(Q/2) 
V2 2 / 




In the case of the triangular lattice this results in a hole energy of 
£if\k*,Q^F) = -5/4 Si with k* = (0,27r/A/3). This is of course only an 
upper limit for the true ground state energy of one single hole; the Trugman 
energy —zi [27] provides the best lower limit known. 
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3 Honeycomb lattice 



Besides the square lattice the honeycomb lattice (see Fig. 9) is one of the few 
examples of bipartite lattice structures in two dimensions. The honeycomb 
lattice is not a Bravais lattice, however, but a triangular lattice with a basis 
of two lattice points. This results in two energy bands with a dispersion 
closely related to the dispersion (2) of the triangular lattice: 



The density of states of the honeycomb lattice is therefore given by 



where Pa{£) stands for the density of states of the triangular lattice discussed 
in section 2. As is shown in Fig. 6 the two bands are separated by a point of 
vanishing density of states at e = 0. The logarithmic van Hove singularities 
correspond to fillings of w = 3/8 and n = 5/8, respectively. 

In the following we are going to investigate the stability of the Nagaoka 
state with respect to a single spin flip on the honeycomb lattice applying 
some of the variational criteria developed in section 2.2. We make use of the 
particle-hole symmetry and consider only the case of less than half filling. Of 
course the tendency towards fully polarized ferromagnetism is expected to 
be much weaker than on the non-bipartite triangular lattice (for more than 
half filling). Compared to the square lattice, which is bipartite too, there 
are nevertheless two important differences. Since the coordination number 
Zh = S is smaller than zd = 4: there exist less loops in the case of the 
honeycomb lattice and the length of the loops is six (see Fig. 9) in contrast 
to a loop length of four on the square lattice. This is important because 
the loops on a lattice are essential for the Nagaoka mechanism as we have 
already stressed in section 1. The other interesting feature is the quasi-gap in 
the density of states at quarter filling (for the Nagaoka state) whereas pn(e) 
shows a logarithmic singularity at this point. 

We have calculated the spin flip energy for the SKA Gutzwiller wave 
function and for its extension by Basile-Elser hopping to nearest neighbour 
sites of the down spin position. As already mentioned in section 2.2 it is 
possible to express all quantities which appear in the spin flip energy (like 
the hole density or the Nagaoka energy) in terms of integrals over the density 
of states and therefore the calculations are completely analogous to those 
for the triangular lattice if one replaces Pa{£) by pH^^)- Fig. 7 shows the 
regions of stability of the Nagaoka state in the phase diagram obtained with 



4(k) =t{M- eA(k)) . 



(17) 




(18) 
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the different variational wave functions. As for the triangular lattice the 
simple Hartree single spin flip ansatz widely overestimates the stabihty of the 
ferromagnetic state. Including the Gutzwiller projector (SKA) we achieve a 
critical hole density of bcv = 0.802 at C/ = oo. The Basile-Elser ansatz for 
nearest neighbours splits the ferromagnetic region into two parts. We find 
an area of possible Nagaoka stability near half filling for large values of U 
(similar to the Basile-Elser instability curve on the square lattice [14]) and a 
stability island around quarter filling, again for large U . Outside of these two 
regions the Nagaoka state is definitely unstable. The critical hole densities at 
infinite V are given by dc:r\ = 0.383 and 8^2 = 0.478 for the instability gap; 
the upper critical hole density is S„3 = 0.662. The stability island extends to 
a minimum critical coupling of Ucr = 16.0 t at ^ = 0.535. As we have already 
observed in Fig. 5 for the triangular lattice the strongest Nagaoka instability 
criterion for small hole densities is again the global instability against phase 
separation between the Nagaoka phase and an antiferromagnetic phase (see 
Fig. 7). 

We have applied the Basile-Elser variational criterion to many lattices 
in different dimensions [14, 28, 29] but only for the honeycomb lattice we 
obtained two separate regions in the phase diagram {0 < 6 < 1) where 
the Nagaoka state is possibly stable against a single spin fiip. What is the 
physical background of this interesting result? First we consider the limit 
of small hole densities. The extension of the Nagaoka theorem by Tasaki [5] 
makes transparent that to achieve a sufficient condition for the existence of a 
fully polarized ferromagnetic ground state in the Nagaoka case (that means 
a half filled system plus one particle at infinite U) each lattice site has to 
be contained in a loop of length three or four. The triangular lattice and 
the square lattice satisfy this criterion whereas the honeycomb lattice does 
not. The possible lack of a Nagaoka theorem suggests that the tendency 
towards ferromagnetism at small hole densities should be even weaker than 
for the square lattice. The exact diagonalization of finite honeycomb chains 
with up to 30 sites by Hirsch [29] shows indeed a degeneracy of the ground 
state energies for different values of the total spin in the small-^ region of 
positive spin fiip energy. As a consequence the Nagaoka state is probably 
not the unique ground state for small hole densities although our variational 
estimates of the spin flip energy are not sufficient to establish the instabihty. 

In the context of our variational calculations the stability of the Nagaoka 
state around quarter filling is an exclusive effect of the quasi-gap in the 
density of states. When 6 approaches 0.5 the Fermi energy drops very quickly 
with increasing hole density thereby enhancing the loss of spin up kinetic 
energy due to the single spin fiip which always contains a term —£p because 
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a spin up particle is removed from the Fermi edge (see section 2.2). This of 
course raises the stability of the Nagaoka state. Figure 8 shows the spin flip 
energy for the SKA Gutzwiller ansatz and for the Basile-Elser ansatz with 
nearest neighbour hopping at C/ = 00, where the effect explained above is 
clearly recognizable. The results of exact diagonalization [29] definitely show 
a maximum total spin of the ground state around 8 = 0.5 thereby suggesting 
that the stabihty island for the Nagaoka state resulting from our calculations 
is real. While we can exclude Nagaoka ferromagnetism in wide parts of the 
phase diagram there is a high probability for a fully polarized ferromagnetic 
ground state for S ~ 0.5 and large values of U. 

4 Kagome lattice 

The so-called kagome lattice (kagome is the Japanese word for basketwork) is 
the line graph of the honeycomb lattice. Generally a line graph is constructed 
by placing a vertex in the middle of each bond of the original lattice and then 
connecting two of the new vertices if the corresponding bonds have a common 
vertex in the original lattice [7]. It follows immediately from this construction 
principle that starting from a lattice with z bonds emanating from each vertex 
the coordination number of the line graph is given by 2 = 2 (2 — 1). It is 
demonstrated in Fig. 9 how the kagome lattice is formed following the rules 
explained above. Since the lattice structure contains triangles it is obvious 
that the kagome lattice is not bipartite in contrast to the honeycomb lattice. 
It can, however, be decomposed into three sublattices each of them having 
again the kagome structure. 

The kagome lattice is not only of academic interest because the kagome 
structure is found in real systems like the anisotropic antiferromagnet 
SrCr8_;sGa4_|_a;Oi9 where the Cr^"*" ions (with spin 3/2) are located in kagome 
planes [16]. Moreover the line graphs of other two-dimensional lattices like 
the square lattice or the triangular lattice show very complicated structures 
due to their high coordination numbers while the kagome lattice has only 
2 = 4. Finally it was the kagome lattice for which Mielke presented the first 
thermodynamically relevant proof of the stability of the Nagaoka state in a 
Hubbard model [7] . The proof is based upon a common feature of all line 
graphs of planar lattices which is the existence of a fiat energy band located 
at the upper band edge for positive t and at the lower band edge for negative 
t, respectively. This situation already makes transparent that especially in 
the case of negative t a very strong tendency towards ferromagnetism is to be 
expected because it is possible to fill particles with the same spin direction 
into the fiat band without any additional cost in kinetic energy compared to 
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a paramagnetic state. 

For t = — 1 the general relation between the densities of states of the 
original lattice (LA) and the line graph (LG) is given by [30] 

Plg{£) = a S{£ + 2) + {1- a) Pla{£ - z + 2) (19) 

where a stands for the weight of the flat band appearing at e = —2. A gap 
between the flat band and the lower edge of the continous spectrum only 
exists if the original lattice is non-bipartite. The upper band edge is given 
by the coordination number z of the line graph. Changing the sign oft means 
replacing e by —e in (19) such that the flat band (and possibly the gap) then 
appear at the upper band edge. 

If |y| = A/" is the number of vertices of the original lattice and \E\ the 
number of edges which is equal to zN/2 the weight of the flat band is given 

by [7] 



\E\ 




\V\ + a 




\E\ 





with (T = 1 if the original lattice is bipartite and otherwise. Making use of 
this result and of equation (18) we obtain the density of states of the kagome 
lattice {t = —1) as 

PK{e) = ^S{e + 2) + ^\s-l\ {{s - if - s) . (21) 

This means that px is composed of a delta peak at £ = — 2 with weight 1/3 
and the two bands of the honeycomb lattice shifted in energy by +1 and 
with an amplitude reduced by a factor of 2/3. Since the honeycomb lattice 
is bipartite no gap appears in this case. 

As already mentioned above, for t = — 1 the fully polarized ferromagnetic 
state obviously is among the possible ground states of the Hubbard model 
for low particle densities (0 < n < 1/3), i. e. in the regime where the Fermi 
energy is in the flat band. In this regime the Nagaoka state is not the 
unique ground state, however, and the degeneracy of the ground state grows 
exponentially with the system size. Mielke [7] used statistical arguments to 
show that at least in part of this density regime the Coulomb energy restricts 
the total spin to be extensive and therefore the system is really ferromagnetic. 

We have applied variational wave fimctions to the kagome lattice in order 
to investigate the instability of the Nagaoka state with respect to a single 
spin flip over the whole density range < n < 2. For our calculations we 
again assumed < n < 1, allowing both signs for the hopping matrix element 
t. First we want to consider the case of more than half filling (t = 1 and 
1 < n < 2) which is mapped to t = —1 and < n < 1 by particle-hole 



20 



transformation such that we recover the situation with the flat band at the 
lower band edge discussed above. As long as all electrons are flUed into the 
flat band (i. e. for low particle densities) the Nagaoka energy is just given 
by e-N = and the local stability of the fully polarized ferromagnet for 
arbitrary U is obvious. The spin flip energy for the SKA ansatz at C/ = oo 
in the regime of the flat band is given by the simple expression 

with the hole density d = \ — n. Since the spin flip energy is also positive 
for smaller hole densities the SKA Gutzwiller ansatz doesn't provide any 
instability of the Nagaoka state at infinite U . Our variational result is of 
course not exact for ^ > 2/3: the true spin flip energy in the regime of the 
flat band is zero. The phase diagram (Fig. 11) shows a strong tendency 
towards fully polarized ferromagnetism also for flnite U . Only for small U 
and small hole densities, i. e. near half flUing, we find a region of guaranteed 
instability of the Nagaoka state. At b = 1/3 the stability of the Nagaoka 
state is slightly enhanced due to the quasi-gap in the density of states taken 
over from the honeycomb lattice (see section 3). The critical coupling Ucr 
goes to zero as the hole density approaches 2/3 which is consistent with the 
stability of the Nagaoka state in the regime of the flat band found by Mielke 
and Tasaki. The results obtained with the SKA Gutzwiller wave function are 
only slightly improved by the Basile-Elser ansatz for nearest neighbours in 
contrast to what happened for the triangular lattice where this ansatz lead to 
a considerable reduction of the Nagaoka stability region. The reason is that 
the additional spin-up hopping processes are much less important if V and 
b are small because the energy cost of double occupancies is moderate and 
the probability for flnding the neighbouring sites of the spin-down position 
already occupied is quite high. 

As explained in section 1 one would not expect to flnd any Nagaoka fer- 
romagnetism in the phase diagram for less than half filling since the kagome 
lattice is not bipartite. With the fiat band now appearing at the upper band 
edge £t = 2 the Nagaoka energy for small hole densities ((5 < 1/3) is given by 
e}4 = —£i 6. As a consequence the loss of spin up kinetic energy compared to 
the Nagaoka state due to a single spin flip is strictly zero because the Fermi 
energy is equal to Sf and the Nagaoka state is therefore locally unstable in the 
whole regime of the fiat band. The spin flip energy for the SKA Gutzwiller 
ansatz (C/ = oo) is calculated to be 

Aeoo = = -3 (5 < (0<(5<l/3). (23) 
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This result holds also for the Basile-Elser ansatz for nearest neighbours be- 
cause the spin-up kinetic energy cannot be further improved. Fig. 10 shows 
the spin flip energy for both variational wave functions at infinite U. The 
quasi-gap in the density of states now appears at ^ = 2/3 and leads again to 
a pronounced enhancement of the loss of spin up kinetic energy compared to 
the Nagaoka state. This effect is strong enough to produce a region of posi- 
tive spin flip energy around 6 = 2/3 even for the Basile-Elser wave function. 
The Nagaoka state is definitely unstable for 6 < 0.672 and for 6 > 0.727 at 
U = oo. With decreasing U the stability area shrinks and finally disappears 
at Ucr = 129.5 t. Up to now there exist only very few examples where the 
variational criteria of SKA and Basile-Elser leave a region of possible stabil- 
ity of the Nagaoka state for n < 1 on non-bipartite lattices. The /cc-lattices 
in dimensions d = 4 and d = 5 show extremely narrow ferromagnetic regions 
for n > 0.99 and large U [29]. We expect, however, that in this limit the Na- 
gaoka state is globally unstable due to the results we presented in section 2.4. 
The situation in the case of the kagome lattice is qualitatively different and 
quite unique because we find a stabihty island well away from half filling for 
n < 1. We cannot exclude, however, that this result might be influenced by 
the fact that we used quite simple variational wave functions where only the 
density of states enters in the calculation of the spin flip energy. A numerical 
check by exact diagonalization is probably difficult because the system size 
has to be quite large in order to reach the flat band situation. 

We have investigated the Nagaoka instability against phase separation 
for the kagome lattice, too. Due to the high ground state degeneracy of 
the classical Heisenberg model on the kagome lattice, it is not at all obvious 
how to choose the magnetic structure for which the spin density calculation 
at half flUing is performed. In order to deal with a simple structure we 
took the special choice of the so called q = phase where the basis cell 
consisting of three spins orientated in a 120° structure is repeated over the 
whole lattice. The phase separation region obtained this way is included in 
the phase diagram of Fig. 11. 

5 Conclusions 

We have investigated the instabifity of the fully polarized ferromagnetic state 
(Nagaoka state) on the triangular, honeycomb and kagome lattices. In com- 
bination with a previous analysis of the square lattice [14] this study gave an 
insight into how the possible stability range of the Nagaoka state depends on 
the lattice structure. In particular, we verifled that more than half-fllled non- 
bipartite lattices generally are very promising candidates for the existence of 
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a Nagaoka region. For less than half filling we were able to exclude a fully 
polarized ferromagnetic ground state on the triangular lattice and reduced 
the possible Nagaoka stability region on the kagome lattice to a small area for 
large U around the filling n = 1/3. We suppose that this behaviour is related 
to the existence of loops of length three in these lattices and the necessarity 
of constructive interference if a particle or hole moves once around a loop. 

Our investigation of the triangular lattice started with a HF approxima- 
tion which gave a survey over the possible magnetic phases. We concluded 
from the missing convexity of the energy as a function of the density for 
more than half filling that the system is unstable with respect to phase sepa- 
ration between the antiferromagnetic 120° structure (favoured at half fiUing) 
and the Nagaoka phase. Evaluating this instability criterion we achieved 
the most pronounced reduction of the possible Nagaoka stability range for 
a particle density slightly above half filling. Below half filling we recognized 
that the 120° structure remains a possible ground state and competes with 
the Nagaoka state. Therefore we suppose that in this case the instability of 
the Nagaoka state might be of global character. We analysed this question 
in more detail presenting an intuitive but strict argument that for arbitrary 
non-bipartite lattices with two or more band maxima at least for small hole 
densities S = 1 — n there are helical states with an energy gain per hole 
exceeding the Nagaoka result by a finite amount. 

We continued our study of the triangular lattice performing a systematic 
variational analysis of the instabihty of the Nagaoka state with respect to 
a single spin flip. This instability turns out to be irrelevant due to other 
instabilities in the framework of HF theory. The inclusion of a Gutzwiller 
projector (SKA), however, allows us to exclude Nagaoka ferromagnetism all- 
together below half filling and leads to a significant improvement of the HF 
estimates for more than half filling. Further improvement is achieved by in- 
cluding majority spin correlation effects into the trial wave function, namely 
hopping processes away from the position of the fiipped spin as first intro- 
duced by Basile and Elser. An efficient calculation of the spin fiip energy for 
this variational ansatz was possible via the extensive use of complete elliptic 
integrals. We showed that even for infinite Coulomb repulsion the Nagaoka 
state is not completely stable for more than half filling, but becomes unstable 
above a critical density = 1.887. Allowing the spin fiip wave functions to 
gain the antiferromagnetic exchange energy we were able to make the critical 
coupling U„ diverge in the limit n \ 1 . The same type of behaviour was 
also achieved by evaluating the physically different approach of localizing 
the flipped spin at a flxed lattice site and solving a single particle scattering 
problem as originally suggested by Richmond and Rickayzen. For this ansatz 
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we obtained, in fact, the best asymptotic resiilt in the limit of n \ 1 for a 
local instability of the Nagaoka state. As mentioned above, however, the 
instabihty is dominated by phase separation in this limit. 

For the bipartite honeycomb lattice we found a Nagaoka instabihty for 
infinite U above a critical hole density S„ = 0.802 already by evaluating 
the SKA Gutzwiller trial wave function. We again improved this ansatz by 
Basile-Elser hopping to nearest neighbour sites and obtained the interesting 
result that the Nagaoka stabihty region splits into two parts, one situated 
near half filling and the other around quarter filling. The upper critical hole 
density is reduced to S„ = 0.662. We traced the probable stability of the 
Nagaoka state around quarter filling - confirmed by finite cluster results - to 
the existence of a quasi-gap in the density of states. The Nagaoka stability 
against a single spin flip near half fllling might be irrelevant since flnite cluster 
results suggest a global instability with respect to a paramagnetic state in 
this region. 

For the kagome lattice we demonstrated that for more than half filling 
there is a strong tendency towards Nagaoka ferromagnetism beyond the flat 
band region 5/3 < n < 2 considered by Mielke. Consistent with the existence 
of the flat band the critical coupling Ucr goes to zero as n approaches 5/3. In 
contrast to the results obtained for other non-bipartite lattices (triangular, 
fee) we found a small ferromagnetic region for very large U even below half 
filling. In the framework of our variational single spin flip states this result 
is again a consequence of the quasi-gap in the density of states occuring at 
n = 1/3 strongly favouring the stability of the Nagaoka state. 

In order to achieve further improvement of the variational single spin flip 
results there exist many possibilities to make the trial wave functions more 
flexible, e.g. by including majority spin correlation effects as it was already 
done by two of us for the square lattice [14]. Furthermore it might be helpful 
to give up the restriction that the removed up spin has a deflnite momentum. 
Investigations of SKA and Basile-Elser spin wave states for multiple chains 
and the square lattice [31] show a significant reduction of the energy due to 
the formation of a bound state between the hole and the flipped spin. 
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Appendix 

The single particle density matrix of the tri- 
angular lattice 

The elements of the single particle density matrix of the Nagaoka state can 
be expressed in the following way: 

h{l) = (co^ct) = f p^^n{£)d£ (A.l) 

where the integers m and n represent coordinates of the vector 1 with respect 
to primitive lattice vectors: 1 = (m + n/2, \^n/2^. The integration extends 
from the Fermi energy to the upper band edge. The functions pm-ni^) are 
generalized densities of states and for the triangular lattice they take the 
following form: 

P™,.(£) = (e*' <?(£ - eA(k)))Bz (A.2) 
= ^cos |^(2m + n)y^ cos (5(£ - £A(k))^ 

The momentum average has to be taken over the first Brillouin zone. For 
m = n = (A.2) obviously describes the usual density of states Pa{£) and 
h{\ = 0) is the hole density 6. The point symmetry of the triangular lattice 
implies 

We therefore need pm,n only for m,n > 0. 

Using Moivre's formula for 171,12 > 0, pm,n can be written as: 

['m+n/2i t ln/2\ j 

p™,™= e E E E(-i)^^''C™r)(;)fe)0T2™+._2„»-2,. (A.4) 



Here ^ stands for 



TuAe) = /cos'" I ^ 1 cos 



^f4^V(^-^A(k))\ (A.5) 



2 / V 2 / / 



BZ 



with u and v being non-negative integers, u + v even. We eliminate the ky- 
integration through the (^-function. Then the t^^v attain the form of elliptic 
integrals in k^. In the following t = — 1 will be assumed. Carrying out the 
substitution x = cos(A;a;) one obtains 

ru,. = -feS^ E (1) (-1 - ^/2)^-^ (A-6) 
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with 



{x + l)-'dx 



The denominator ^{x) is the positive root of 

e{x) = ii-x'){s+E-{x-i-^Er] 



(A.7) 



(A.8) 



Concerning the integration limits in (A.7) the two cases of £ > — 2 or e < — 2 
have to be considered separately. For values of e less than —2 the integration 
hmits are given by 



X^ = 1 + -£ — \/3 + £ 



X+ = 1 + -£ + V3 + £ . 



(A.9) 



If £ > — 2 the upper limit is 1. 

It follows from the construction of the t^^v (A. 6) that the index s of J, can 
be any integer. For negative indices the binomial theorem allows to replace 
the (s > 0) by other elliptic integrals 



^-. = E(:)^^(^>o) 

r=0 



with 



J X- 



H^+'i) x^'dx 



(r > 0) . 



(A.IO) 



(A.ll) 



By integrating -£:{x + l)~''^(a;) and -£:X^^{x), respectively, from a;_ to 
min(a;+, 1), which are both zeros of the function ^(a:), one obtains the follow- 
ing two recursion formulae [32] : 



J, = 



(2s - l)bi 
+2(2 - 5)64^.-3 



2sboJs+i + 2(1 - s)b2Js-i + (3 - 2s)63J.-2 

s > 1) , (A.12) 
2s - 3)a3J,_i + 2{s - 2)a2ls-2 + (2s - 5)aiJ,_3 

(A.13) 



2(1 - s)a^ 

+2(5 - 3)ao/,_4l (5 > 3) . 



The parameters and hi are the coefficients obtained if one writes i^{x) as 
a polynomial in x and a; + 1, respectively: 



e{x) = j2a^x' = j2h^{x^ly . 

i=0 i=0 



(A.14) 
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Note that for 5 = 3 the term containing /_i in (A. 13) vanishes. The term 
containing J^+i in (A. 12) vanishes for all s since feo = 0. Due to these 
two peculiarities the recursion formulae (A. 12) and (A. 13) allow to express 
all Is and by Jq, h and I2. These three functions are given by linear 
combinations of complete elliptic integrals of the first and second kind K[m] 
and E[m], respectively (see [32]): 

-z, E [zllz^^ + ^yiA + 2 + K [zH zl] 

(max(a;+, 1) — a;_)(min(a;+, 1) + 1) , (A. 15) 

(max(a;+, 1) + l)(min(a;+, 1) — a;_) , 
J i(l + £/2)2 for £ < -2 , 

j 2 + e/2-^/3T7 for e > -2 . 

Using the above formulae we see that the generalized densities of states 
are all given by linear combinations of complete elliptic integrals of the first 
and second kind. For the calculation of the elements of the single particle 
density matrix only one numerical integration has to be performed. 

Equations (A. 4) and (A. 6) show that the density of states pa (see Fig. 1) 
is proportional to Jq and therefore can be expressed by a complete elliptic 
integral of the first kind: 

PA = — i^UM] • (A.16) 

All the other p^^^ are most conveniently obtained by implementing the above 
formulae using an algebraic manipulation program. 



h = 
h = 
with Zq = 
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Figure captions 



Figure 1 Density of states of tlie triangular lattice for i = 1. 

Figure 2 HF phase diagram of the triangular lattice for t > 0. The abbrevi- 
ation AF marks the stability region of the antiferromagnetic 120° structure, 
N stands for the colhnear Neel structure, 5*1 for the phase separation re- 
gion AF-Nagaoka, whereas ^2 corresponds to the phase separation region 
AF-Helical. Hatching marks the regions where we find helical states with 
continuously varying Q vectors. 

Figure 3 Brillouin zone of the triangular lattice and illustration of the path 
followed by the favoured Q vectors along the para-helical phase boundary 
with increasing density < n < 1.5. 

Figure 4 Phase diagram for the triangular lattice: Nagaoka instability 
lines obtained with the Hartree single spin ilip ansatz (dotted), the SKA 
Gutzwiller ansatz (long dashed), the Basile-Elser ansatz with nearest neigh- 
bour hopping (short dashed) and the best Basile-Elser like wave function 
including the exchange term (14) (full line), t = 1. Hatching marks the 
region where the Nagaoka state might be a ground state. 

Figure 5 Phase diagram showing the instability of the Nagaoka state on the 
triangular lattice in the limit n \ 1: the instability lines for the Richmond- 
Rickayzen ansatz (full line) and for the single spin flip wave function including 
the exchange term (14) (dotted) are compared to the HF criterion for the 
instability against infinitesimal twist (long dashed) and the global instability 
against phase separation (short dashed). 

Figure 6 Density of states of the honeycomb lattice for \t\ = 1. 

Figure 7 Phase diagram for the honeycomb lattice: Nagaoka instability 
lines obtained with the Hartree single spin flip ansatz (dotted), the SKA 
Gutzwiller ansatz (short dashed) and the Basile-Elser ansatz with nearest 
neighbour hopping (full line). The global instability against phase separation 
(long dashed) is also shown. Hatching marks the regions where the Nagaoka 
state might be a ground state. The phase diagram for more than half fllling 
is obtained by particle-hole transformation which means replacing 6 with n. 
Ubr = 12.59|i| 

Figure 8 Spin flip energy Agoo on the honeycomb lattice at inflnite U for 
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the SKA Gutzwiller ansatz (dashed) and the Basile-Elser ansatz with nearest 
neighbour hopping (full line). The quasi-gap in the density of states (Fig. 6) 
corresponds to ^ = 0.5. 

Figure 9 Construction of the kagome lattice as the line graph of the honey- 
comb lattice. 

Figure 10 Spin flip energy Acoo on the kagome lattice at infinite U for the 
SKA Gutzwiller ansatz (dashed) and the Basile-Elser ansatz with nearest 
neighbour hopping (full line), for less than half filling {n < 1). The regime 
< (5 < 1/3 corresponds to the flat band at the upper band edge. 

Figure 11 Phase diagram for the kagome lattice: Nagaoka instability hues 
obtained with the Hartree single spin flip ansatz (dotted), the SKA Gutzwiller 
ansatz (short dashed) and the Basile-Elser ansatz with nearest neighbour 
hopping (full line). The global instability against phase separation (long 
dashed) is also shown. The regimes 2/3 < n < 1 and 5/3 < n < 2 cor- 
respond to the flat bands. Hatching marks the regions where the Nagaoka 
state might be a ground state. Ubr = 13.80|t| 
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Figure 1: Density of states of the triangular lattice for t = 1. 
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Figure 2: HF phase diagram of the triangular lattice for t > 0. The abbrevi- 
ation AF marks the stability region of the antiferromagnetic 120° structure, 
N stands for the coUinear Neel structure, 5*1 for the phase separation re- 
gion AF-Nagaoka, whereas ^2 corresponds to the phase separation region 
AF-Helical. 
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Figure 3: Brillouin zone of the triangular lattice and illustration of the path 
followed by the favoured Q vectors along the para-helical phase boundary 
with increasing density < n < 1.5. 
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Figure 4: Phase diagram for the triangular lattice: Nagaoka instability 
lines obtained with the Hartree single spin flip ansatz (dotted), the SKA 
Gutzwiller ansatz (long dashed), the Basile-Elser ansatz with nearest neigh- 
bour hopping (short dashed) and the best Basile-Elser hke wave function 
including the exchange term (14) (full line), t = 1. 
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Figure 5: Phase diagram showing the instability of the Nagaoka state on the 
triangular lattice in the limit n \ 1: the instability lines for the Richmond- 
Rickayzen ansatz (full line) and for the single spin flip wave function including 
the exchange term (14) (dotted) are compared to the HF criterion for the 
instability against infinitesimal twist (long dashed) and the global instability 
against phase separation (short dashed). 



36 




-3-2-10 1 2 3 



energy 

Figure 6: Density of states of the honeycomb lattice for \t\ = 1. 
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Figure 7: Phase diagram for the honeycomb lattice: Nagaoka instability 
lines obtained with the Hartree single spin flip ansatz (dotted), the SKA 
Gutzwiller ansatz (short dashed) and the Basile-Elser ansatz with nearest 
neighbour hopping (full line). The global instability against phase separation 
(long dashed) is also shown. The phase diagram for more than half fllling 
is obtained by particle-hole transformation which means replacing 6 with n. 
Ubr = I2.59\t\ 
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Figure 8: Spin flip energy Aboo on the honeycomb lattice at infinite U for 
the SKA Gutzwiller ansatz (dashed) and the Basile-Elser ansatz with nearest 
neighbour hopping (full line). The quasi-gap in the density of states (Fig. 6) 
corresponds to ^ = 0.5. 
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Figure 10: Spin flip energy Acoo on the kagome lattice at inflnite U for the 
SKA Gutzwiller ansatz (dashed) and the Basile-Elser ansatz with nearest 
neighbour hopping (full line), for less than half filling (n < 1). The regime 
< (5 < 1/3 corresponds to the flat band at the upper band edge. 
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Figure 11: Phase diagram for the kagome lattice: Nagaoka instability 
lines obtained with the Hartree single spin flip ansatz (dotted), the SKA 
Gutzwiller ansatz (short dashed) and the Basile-Elser ansatz with nearest 
neighbour hopping (full line). The global instability against phase separa- 
tion (long dashed) is also shown. The regimes 2/3 < w < 1 and 5/3 < w < 2 
correspond to the flat bands. Ubr = 13.80|i| 
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